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High-resolution optical imaging of the core of the globular 
cluster M15 with FastCam 
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ABSTRACT 

We present high-resolution /-band imaging of the core of the globular cluster M15 
obtained at the 2.5 m Nordic Optical Telescope with FastCam, a low readout noise 
L3CCD based instrument. Short exposure times (30 ms) were used to record 200000 
images (512 x 512 pixels each) over a period of 2 hours 43 min. The lucky imag- 
ing technique was then applied to generate a final image of the cluster centre with 
FWHM~ 0".l and 13" x 13" FoV. We obtained a catalogue of objects in this re- 
gion with a limiting magnitude of I = 19.5. /-band photometry and astrometry are 
reported for 1181 stars. This is the deepest /-band observation of the M15 core at 
this spatial resolution. Simulations show that crowding is limiting the completeness of 
the catalogue. At shorter wavelengths, a similar number of objects has been reported 
using HST/WFPC observations of the same field. The cross-match with the available 
HST catalogues allowed us to produce colour-magnitude diagrams where we identify 
new Blue Straggler star candidates and previously known stars of this class. 

Key words: instrumentation: high angular resolution - techniques: photometric - 
binaries: close - stars: blue stragglers - globular clusters: individual: M15 - Hertzsprun- 
gRussell and colourmagnitude diagrams 
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1 INTRODUCTION 

The cores of globular clusters contain a very dense stellar 
population reaching up to 10 6 stars pc -3 . They are test 
laboratories for studying the dynamics in dense environ- 
ments and provide crucial information on the evolution of 
the cluster. Encounters and direct collisions between core 
cluster stars are highly probable leading to a variety of ex- 
otic species of stellar objects such as close interacting bina- 
ries, pulsars, x-ray sources and Blue Stragglers (BS). The 
globular cluster cores are also potential hosts of intermedi- 
ate mass black holes (|McNamara et all (|2003l )h which may 
be revealed by the dynamics of the stars in their very inner 
regions. Resolving the core stellar population is required to 
establish the existence and determine the mass of any cen- 
tral compact object. Studies of the stellar populations in 
the globular cluster most inner regions have been conducted 
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in the optical usin g the Hubble Space Telesco pe (see e.g. 
lYannv et all l| 19941 ). iGuhathakurta et al.1 1| 19961) ) and in the 
near infrared using Adaptive Optics systems in g round based 
telescopes (see e.g. iDavidge fc Courteaul |l99^)). 

M15 is one of the oldest and most massive globular 
clusters with an unus ually compact core and a high cen- 
tral velocity dispersion (|McNamara et all (|2003l )h This clus- 
ter is particularly suitable to investigate the presence of 
a massive central compact object. High-resolution obser- 
vations of t h e M1 5 core have already been obtained by 
lYannv et all (|l994h with the HST/PC camera with the fil- 
ters F336W, F555W an d F785LP (UVI bands) and by 
iFerraro fc Parescd ^199$ ) using F140W, F220W and F342W 
filters and the HST/FOC camera. These early data were 
taken with aberrated optics and it was difficult to per- 
form accurate photometry of the core stars. After post- 
repair HST took images of the core with the HST/WFPC2 
with the filters F336W, F4 39W and F555W (UBV bands) 
|Guhathakurta et al.l |l996h b and more data in the V band 
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were reported by McN amara et all l|2003h . High resolution 
optical images using adaptive optics in the Canada- France- 
Hawa ii Telescope have be en r eported by iGebhardt et"al] 
(|2000h . iPiotto et all (|2002h and Ivan der Marel et all (|2002l ) 
also performed independen t photometry of the data taken by 
iGuhathakurta et alj (119961 ) and in particular, the catalogue 
by Piotto et alj (I2002T) has been used to study Blue Strag- 
gler stars bv lMoretti et all (120081 ). The inner region of M15 
has been studied more recently using deep FUV F140LP and 
NUV F220W i mage s obtained with the HST/ACS camera 
(jPieball et al.l l|2007l )) which show the presence of horizontal 
branch stars, blue stragglers and white dwarfs. 

The development of fast readout very low noise L3CCD 
detectors has opened the possibility of obtaining high- 
resolution images of the cores of globular clusters using 
ground based tel escopes and the lucky imaging technique 
l|Law et all (|2006l )). Here we present lucky imaging observa- 
tions in the I-band of th e core of the globul ar cluster Ml 5 
obtained with FastCam (jOscoz et al] (|2008l )) at the 2.5m 
Nordic Optical Telescope. We report /-band photometry for 
1181 sources in a radius of ~ 6.5 arcsec from the centre 
of M15 and compare with available stellar catalogues from 
HST. We find that the lucky imaging technique provides re- 
liable photometric measurements for stars in crowded fields. 
We take advantage of this technique to study known variable 
stars and to search for new objects in the core. The same 
instrument was used recently to obt ain high contrast opt i- 
cal imaging of substellar companions (jLabadie et all (|201ll )) 
and combined with adaptive optics t o produce high precisio n 
astrometry of a brown dwarf binary (|Femema et al] l|201ll) ). 



region for each image. We selected the 70 images with the 
highest values of the brightest pixel (i.e. 7 % of the total). 
We determined the position of the brightest pixels and ap- 
plied a shift-and-add algorithm to these selected images and 
generated one final image per cube. Our shift-and-add algo- 
rithm employs integral pixel shifts (sub-pixel shifts provide 
a slight improvement of order 5 % of the FWHM of the final 
image with respect the case of integer shifts). In this way, 
we obtained 200 images with an effective integration time of 
2.1 s each. This procedure was carried out for different per- 
cent criteria and found that the best compromise between 
sensitivity and spatial resolution was achieved by selecting a 
fraction of images in the range 5-10%, therefore we adopted 
the 7% criterion. 

As a second step we took a reference image from these 
200 images and defined a new circular region of 45 pixels of 
diameter centred in the previous bright star. Those images 
where the selected bright star was out of this region were 
ruled out in order to avoid too large shifts for the final com- 
bined image. Only two images were ruled out of 200. Then 
we applied the shift-and-add algorithm again and obtained a 
final image of 512 x 512 pixels. Only in the central 422 x 422 
pixels there is an effective integration time of 415.8 s, due 
to the maximum shift applied to the images (512 — 2 x 45 
pixels) . This central field is shown on the left panel of figure 
[T]and we will focus on this region to carry out the photomet- 
ric work. On the right panel we display observations of the 
same field with HST/WFPC2 ( centra l chip) in the F555W 
(V) filter (|Guhathakurta et al] (| 19961 )). 



2 OBSERVATIONS AND DATA REDUCTION 

Observations were obtained on the night of 2008 July 25 
with the FastCam instrument on the 2.5-m Nordic Optical 
Telescope using the /-band filter (820 nm) which matches 
the Johnson-Bessell system. FastCam was equipped with 
a 512 x 512 pixels L3CCD from Andor Technology. The 
optics was designed to provide a spatial scale of about 30 
mas/pixel which was later determined accurately via com- 
parison with available HST astrometry for M15 stars and 
using vi sual binaries f rom t he Washington double star cat- 
alogue dMason et al] (|200ll )). We pointed the telescope to 
the centre of M15 (R.A.= 21 h 29 m 58 s .3 and Dec= 12°10'l") 
and collected 200 series ("cubes") of 1000 images each with 
individual exposures of 30 ms. These 200000 images were 
obtained during 2 hours 43 min of real observing time. The 
average seeing during this observing period was ~ 0.65 arc- 
sec as measured in the final /-band image (combining all of 
the I-band images without applying the shift-and-add algo- 
rithm) . 

We made the raw data reduction in two steps. First, we 
selected the best 70 images of each cube (7% of the total) as 
follows. We defined a circular region on the first image of a 
cube with a 45 pixels diamete r ( T'.3), centred in the brig ht- 
est star (ID 6290 in reference Ivan der Marel etafl (|2002n . it 
is marked with a yellow circle on the left panel of figure [TJ 
near the centre of the field. This star was always the bright- 
est available within this circular region. We checked that 
in all the images the reference star was present within this 
region. Then we identified the brightest pixel within this 



3 ANALYSIS 

In Fig. [2] we show the average of the light profile of 14 rel- 
atively isolated stars in our image. Each stellar profile is 
normalised and a PSF template is built with a shift-and- 
add algorithm at sub-pixel scales, each pixel of the PSF 
template is shown with a point in the figure. We note the 
presence of a remarkable narrow core and an extended halo. 
This is the characteristic profile of stars in images obtained 
using the lucky imag i ng tec hnique as described in previous 
work (Bal dwin et all (|2008l )). The FWHM of the average 
normalised PSF in our final image is 0".l very close to the 
diffraction limit of 0".084 of the telescope at the wavelength 
of the observations. The normalised light profiles for two 
stars are shown in figure [3] the profile marked with red cir- 
cles belongs to a star at a distance of 2". 2 from the the star 
that was used for image alignment and the profiles marked 
with crosses (x) belongs to a star at a distance of 5".l. There 
are slight differences between the two which may indicate a 
modest amount of anisoplanicity and that the dominant tur- 
bulent layer at the time of the observations was probably at 
very low altitude. It should be pointed out that the stars 
are not perfectly round. We do not have an obvious expla- 
nation for this, in particular it could be related to guiding 
issues which may be at play. We also note the relatively 
broad azimutal distribution of fluxes in the halo of the two 
stars which is caused by the presence of very close weaker 
contaminant stars. 
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Figure 1. FastCam image of the central region of M15 (left panel) obtained using the lucky imaging technique, the yellow circle shows 
the location of the star that was used for image alignment. The scale plate is 0".031 pixel -1 , and the field of view is 13". 16 X 13". 16. The 
same image after convolution with a "mexican hat" kernel to improve the detectability of point sources (central panel). HST/WFPC2 
(central chip, from Guhathakurta et al. (1996)) in the F555W filter (right panel). North is up and east to the left for all the images. 
They are displayed using a square-root stretch. 
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Figure 2. The average of the light profile of the 14 relatively 
isolated stars in our image. Here the light profile of each star is 
normalised and a PSF template is formed with a shift-and-add 
algorithm at sub-pixel scales, each pixel of the PSF template is 
shown with a point in the figure. 
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Figure 3. The normalised light profiles for two stars: the profile 
marked with red circles belongs to a star at a distance of 2". 2 
from the star that was used for image alignment and the profiles 
marked with crosses (x) belongs to a star at a distance of 5".l. 
Each point in the figures represents a pixel. 



3.1 Detection of stars 

In order to carry out the detection of objects in our final im- 
age we applied a derivative filter to the whole image, we used 
the Laplacian of Gaussian filter, often called " mexican hat" , 
which works efficiently in dense stellar fields |Debrav et al.l 

Gil)), 

where <r m controls the width of the Gaussian and x, y de- 
termine the coordinates of the pixel. In the image convolved 



with the "mexican hat" (with <r m = 1 pixel), the peaks are 
enhanced and the saddle points approach the background, 
so the detectability of the sources is definitely improved. 

The local maxima are detected on the convolved im- 
age using a threshold, set at 3 times the standard deviation 
of the background determined in the less crowded part of 
the image. In order to measure this background we divided 
the image into squares of 24 x 24 pixels and measured the 
standard deviation for each. The lowest of these values was 
adopted as the standard deviation (a) of the less crowded 
part of the image. We note that no region was completely 
free of sources. We determined the local maxima by com- 
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parison of each pixel with its 8 neighbours. In this way we 
generated a list of potential objects with their positions de- 
termined as the centroid of a box of 3 x 3 pixels centred 
around each maximum. On the central panel of figure [1] we 
can see the resulting convolved image. 1682 objects were 
detected in this image of 422 x 422 pixels (13.16 x 13.16 
arcsec 2 ). Only a very small fraction are false detections. 



3.2 Cross-match with HST/WFPC2 data 

The spatial scale of the image was obtained from a 
cross-match with th e avail able HST/WFPC2 catalogue by 
van der Marel et all j2002T ) as described below. We first 
made a linear transformation of the catalogue into our list 
and identified 762 common pairs with V < 19 within a tol- 
erance of 2 pixels of FastCam. We improved the number of 
matches among the two datasets using a second order trans- 
formation from the x-y plane in our image to the AR.A.- 
ADec. plane of the catalogue. The transformation equations 
modelled a zero point difference, a scale change, a rotation, 
a tilt, and second order distortion corrections. These trans- 
formations include 12 unknown coefficients that are solved 
using a least squares method 

AR.A. = a + bAx + cAy + d{Axf + eAxAy + /(Ay) 2 , (2) 

Adecl. = a' + b'Ax + c'Ay+d'(Ax) 2 + e'AxAy+f'(Ay) 2 , (3) 

with AR.A. = R.A. - R.A. cen t, Adecl. = decl. - decl. ccnt , Ax 
= x - x ccnt and Ay = y - y ccn t, where R.A. cont , decl. ccnt , x ccnt 
and y cen t are taken from the object nearest to the centre in 
our image. 

Using the initial list of 762 matched objects with V < 19 
we found the transformation coefficients using least squares. 
These coefficients were used to compute the positions of our 
objects in the AR.A.-ADec. plane. For each object in our 
list, we determined distances to the objects in the HST cat- 
alogue, the minimum of these distances determines the cor- 
relation length for the object. We assume that a FastCam 
object can be matched with an HST object when the cor- 
relation length is smaller than 0.046 arcsec (1.5 pixels). We 
removed from the original list used for the transformation 
objects with correlation length larger than 0.046 arcsec and 
V < 19 and performed an iterative process until we ob- 
tained a collection of FastCam objects which all have a cor- 
relation length smaller than such value. The final list of ob- 
jects used to calculate the coefficients of the transformations 
have 770 objects with V < 19 distributed uniformly over the 
whole image. Then we used these transformations to cross- 
match the 1682 objects fo und by FastCam wit h obje cts in 
the catalogue of reference Ivan der Marel et all (|2002h and 
finally selected those with a correlation length smaller than 
0.046 arcsec, whi ch corresponds to 1 pixel in t he central chip 
of HST/WFPC2 Ivan der Marel et all \200j ). We obtained 
1481 matched stars and 201 unmatched, we note that these 
unmatched stars include some potential false detections in 
the FastCam image. 

In the catalogue given by Ivan der Marel et al.l 
J2002I) and obtain e d fro m the images presented in 
iGuhathakurta et alj (| 19961 ). there are 3015 objects in the 
field observed by FastCam. These authors used a combined 
F336W, F439W and F555W image (UBV bands) to detect 
stars and added manually faint stars with peak brightness 



below the detection threshold of the DAOPHOT/FIND 
routine. From the same se t of observa t ions but using 
only the B and V images IPiotto et all l|2002l ) detected 
2221 objects in this field. We can increase the number of 
detections in our image by decreasing the threshold used in 
subsection 3.1 to detect local maxima. For a 2a threshold 
we obtain 1911 ob j ects m atched with the catalogue by 
van der Marel etafl <|2002T i and 384 unmatched; and at la, 
we obtain 2275 matched and 1452 unmatched objects. In 
this work we restrict the study to the 1682 objects detected 
at 3a level. We removed the false detections among the 
201 unmatched objects by visual inspection and concluded 
that 46 out of these 201 were true stars, which included 
both new stars non reported by Ivan der Marel et all (|2002T ) 
and stars with too large uncertainty in the measured 
positions to perform a reliable correlation. For 41 of these 
stars the difference between the centroid position and the 
location of the maximum signal in the convolved image was 
more than 1 pixel. This difference is mainly caused by the 
presence of a brighter contaminating star in the vicinity. 
There are 5 unmatched stars. Finally, we verified that the 
transformation coefficients associated with tilt and second 
order distortion corrections were negligible, and therefore 
using a linear transformation was sufficient to determine 
a plate scale of 31.18 ± 0.03 mas/pixel for the FastCam 
image. 



3.3 Photometry 

In order t o perform photomet r y we use the technique 
adopted bv IGuhathakurta et al.l (1996) which incorporates 
a set of standard DAOPHOT routines. Firstly we construct 
an empirical PSF by iteration. In each iterative step we 
construct a PSF template of radius 22 pixels by averaging 
14 bright and relatively isolated stars. Then we remove the 
neighbours of these stars using the PSF template from the 
previous iteration and calculate again the PSF template of 
these stars. In each iteration the quality of the PSF is im- 
proved because the neighbours are removed with increasing 
precision, convergence is reached after five iterations. From 
all the template functions tested the penny2 function which 
vary linearly with the image position provides the best ap- 
proximation to the actual PSF. The penny2 function is a 
Gaussian core with a Lorentzian profile, in which both com- 
ponents can be tilted. Then the final PSF template is fitted 
to all stars detected on the FastCam image with the ALL- 
STAR routine in DAOPHOT. Due to the difficulty to find 
relatively isolated stars on the image, the PSF subtraction 
does not work perfectly and the stars are not subtracted 
completely. This is because the PSF used for subtraction of 
the flux of neighbour stars is determined locally using the 
brightest, nearest stars, but sometimes there are no bright 
stars which can help to define the local PSF with very good 
S/N and therefore the correction may not be perfect. This 
prevents to carry out direct PSF photometry in our image. 
Therefore, in order to obtain accurate photometry we fol- 
lowed the hybrid method described bv IGuhathakurta et al.l 
( 1996) which combines direct aperture photometry and PSF 
fitting. This method works as follows: we take a star from 
our list of objects and subtract the PSF to each neighbour- 
ing star within a 44 pixel radius (two times the radius of 
the region used to calculate the PSF). Then aperture pho- 
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Figure 4. Calibrated FastCam magnitude (Ifc) vs - Yanny et 
al.(1994) magnitude (ly) for the 269 stars in common. Dashed 
lines are the ±2tr dispersion, stars within the region defined by 
dashed lines are chosen for the zeropint calibration. 



tometry is obtained for the star using a circular aperture of 
radius 2.5 pixels (0,078 arcsec) with a local sky background 
measured in a surrounding annulus. This procedure was re- 
peated for each detected target and obtained instrumental 
magnitudes, I/ c . 



3.4 Photometric calibration 

The instrumental Ij c magnitudes were converted to the 
Johnson I standard system using stars in the catalogue by 
lYannv et al.l ljl994l ). These authors had transformed their 
F785LP i nstru mental magnitudes to Johnson / system 
jjohnsonl (l966)) and we directly adopted their / magni- 
tudes. We correlated with this catalogue using the algorithm 
described in subsection 13.21 and obtained 269 stars in com- 
mon. The zeropoint of the calibration was obtained via an 
iterative process using only 219 stars in the magnitude range 
12-18. The transformation was found to be linear and no 
colour term was found to be significant (figure 3]). 

The number of FastCam detected stars in 0.25 mag bins 
as a function of I magnitude is represented in figure [5] We 
can see that this number increases until the 19-19.25 mag 
bin, and reach stars as faint as I m 21.5. We removed 82 stars 
from our catalogue which had false photometry due to their 
proximity to much brighter stars. Our catalogue provides 
/-band photometry of 1181 stars with / ^ 19.5, they are 
ordered from low to high I magnitude and it is presented in 
its entirety in the electronic edition of MNRAS. For reference 
the first 10 entries are listed in table [T] 



3.5 Photometric accuracy 

In order to investigate the statistical errors of the photome- 
try in the final FastCam image, we bin the series of 200 im- 
ages into combinations of 40 images each and perform the 
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Figure 5. Number of stars found with FastCam in 0.25 mag bins 
as a function of the I magnitude. 



Table 1. FastCam /-band catalogue of the M15 core 



ID 


AR.A. 


Adec. 


I 


V 


IDvm 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


1 


4.387 


-0.053 


12.03 


13.369 


4113 


2 


2.101 


-2.216 


12.13 


13.696 


5469 


3 


7.188 


-6.430 


12.17 


13.605 


5166 


4 


-4.717 


2.919 


12.21 


13.352 


6041 


5 


-0.110 


-2.527 


12.27 


13.743 


6290 


6 


-2.850 


-1.718 


12.69 


13.977 


6947 


7 


-5.211 


-4.359 


12.72 


13.900 


8408 


8 


5.711 


-0.607 


12.74 


14.072 


3891 


9 


-4.681 


-6.160 


12.77 


13.870 


8777 


10 


4.732 


-6.416 


12.85 


13.952 


5956 



Notes. - Table 1 is presented in its entirety in the electronic edition 
of MNRAS. The first 10 entries of the /-band FastCam photomet- 
ric catalogue described in the text are shown here for guidance 
regarding its form and content. The full catalogue contains 1181 
stars. Col. (1) is the ID number of the star. Cols. (2) and (3) 
give coordinates (R.A., Dec.) for each star, measured in arcsec- 
onds with respect to the M15 cluster star AC 211 (which ac- 
cording to Kulkarni et al. (1990) has R.A.=21 h 29 m 58 s .310 and 
Dec= 12°10'02".85)). Col. (4) gives the / magnitude obtained 
with FastCam and cols. (5) and (6) the V magnitude and the 
corresponding ID from the van der Marel et al. (2002) catalogue. 



same photometric analysis described above. We therefore 
obtain the photometry for each star in our catalogue which 
is bright enough to be detected in each of the 5 images. Fig- 
ure[6]shows the standard deviation of the 5 measurements as 
a function of stellar magnitude. In the magnitude range 12- 
15 the statistical errors are below 0.05 mag. In the plot we 
mark stars that are known to be variable (see section 4). In 
the final combined image we expect these statistical errors 
to decrease significantly and the final error to be dominated 
by the systematic effects resulting from the variation of the 
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Figure 6. Dispersion of the individual photometry of each star 
versus mean I magnitude, for the 5 photometric measurements of 
combinations of 40 images as explained in the text. Variable stars 
from reference Dieball et al. (2007) are denoted by green circles. 
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Figure 7. Photometric errors vs. / magnitude as inferred from 
simulations. The top panel displays the difference between the 
input simulated magnitude and the output measured magnitude 
for each star as a function of magnitude. On the bottom panel 
the mean value and the dispersion for each 0.5 mag bin (each 
simulated image) are shown. 



PSF with position in the image and the residual contamina- 
tion of poorly subtracted neighbouring stars. 

In order to study the impact of systematic errors we 
simulate stars using the previous PSF template obtained 
with the penny2 function which vary linearly with the posi- 
tion in the image. The stars are simulated in different posi- 
tions on a 422x422 pixels image with initial zero background 
using standard DAOPHOT routines. We perform aperture 
photometry using a circular aperture of radius 2.5 pixels and 
find that the average difference between the measured and 
simulated magnitudes is less than 0.015 mag independent 
of the simulated stars. Such differences are smaller than the 
statistical error and could be caused by a non-perfect ade- 
quacy of the PSF template due to spatial variability. 
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Figure 8. Detectability of stars in our final image vs / magnitude 
as inferred from the simulations. 100 stars were simulated at each 
0.5 mag bin. The dashed line is the 100 % detectability curve. 




input 



Figure 9. Position errors (in milliarcseconds) versus / magnitude 
of stars as recovered from the simulations. The top panel displays 
the difference in position between the input and the recovered 
value from the simulated images. In the bottom panel we show 
mean values and dispersions for stars in each 0.5 mag bin. 



Errors associated to crowding were investigated simu- 
lating stars of a given magnitude with a spatially constant 
PSF template. We explored the magnitude range 1=12 to 20 
mag simulating 100 stars for each step of 0.05 mag. In total, 
16 images were generated each containing simulated stars of 
a given magnitude distributed uniformly on the final Fast- 
Cam image with a minimum imposed separation of 4 pixels 
from any two stars (85% of the stars detected by FastCam 
have their nearest neighbours beyond 4 pixels and our de- 
tection algorithm only consider stars with separation larger 
than 3 pixels). Stellar magnitudes were then measured in 
theses images using the hybrid method described above. In 
the top panel of figure[7]we see the difference between the in- 
put and the output measured magnitudes for the simulated 
stars as a function of input magnitude. At the bottom panel 
of this figure we show the mean value and the dispersion for 
each magnitude bin. We conclude that there is no obvious 
bias in our photometric measurements associated to crowd- 
ing. For the brightest objects the contribution of crowding 
to the photometric error appears to be less than 0.04 mag 
while for objects of magnitude I = 18 is of order 0.2 mag. 
Therefore, in our final image errors associated to crowding 
are similar to the statistical errors for / < 18. 

We also give in figure [8] the number of stars detected in 
the simulated images for each magnitude bin. For instance, 
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Figure 10. Colour-magnitude diagram / vs. U — V. The /-band 
magnitude is from FastCam, U and V magnitudes are from van 
der Marel et al 2002. There are 1312 objects correlated between 
both catalogues with / < 20. 

in the 18.5-19 mag bin more than the 80% of the objects 
originally included in the simulation were detected by the 
algorithm described in subsection 13.11 and for the 19-19.5 
mag bin more than 50% were still detected. To our knowl- 
edge this is the deepest and more complete high spatial res- 
olution /-band observation of the M15 core reported so far. 

We also used the previous set of simulations to inves- 
tigate any possible bias in our astrometry measurements. 
The difference between the input and recovered position for 
each artificial star in the simulated images is plotted at the 
top panel of figure [5] in milliarcseconds. The positions are 
measured by the method described in subsection 13.11 The 
mean value and the dispersion for each bin are shown at the 
bottom panel of figure [5] 

4 STELLAR POPULATIONS IN THE CORE 
OF M15 

Figure [TOl displays the colour-magnitude diagram (CMD) / 
vs. U — V for the 1312 st a rs cro ss-matched with the catalogue 
of Ivan der Marel et all (|2002h and with I < 20. We note 
the main-sequence turnoff at I « 19. The CMD V vs. U — 
I shown in figure [TT] offers the advantage of a long colour 
baseline resulting in a more clear separation of the various 
types of stars. These are denoted in the figure with different 
symbols: blue stragglers (□, BS), bright red giant branch 
(V, Bright RGB), (red) horizontal branch (x, HB) and blue 
horizontal branch (+, BHB). The colour criteria used to 
as sign the vario u s stell ar types have been adopted according 
to I Yanny et all j 19941 ) and are not strict. Number counts 
for the various types of stars found in the FastCam final 
full i mage are giv e n in c olumn 2 of table [21 For comparison 
with lYannv et ail |l994t ) we also list in columns 3 and 4 the 
number counts for r < 5.6 arcsec obtained in this work and 
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U-I 

Figure 11. Colour-magnitude diagram V vs. U — I for the M15 
core stars, /-band magnitudes are from FastCam (this work), U 
and V magnitudes are from HST (van der Marel et al. (2002)). 
Blue stragglers (BS), bright red giant branch (RGB), (red) hor- 
izontal branch (HB) and blue horizontal branch (BHB). Blue 
straggler stars given in previous references are marked with filled 
diamonds. Variable stars are taken from Dieball et al. (2007). AC 
211 is the variable star in the BHB zone. 



Table 2. Counts of various stellar types. 



Stellar type 


Counts 


Counts 


Counts 






(r < 5". 6) 


(r < 5".6) 








(Yanny et al. (1994)) 


Bright RGB 


40 


26 


25 


HB 


29 


22 


18 


BS 


10 


7 


7 


BHB 


20 


12 


12 



by these authors. We find no significant differences between 
the two. The CMD I vs. V — / is also shown in figure 1121 
Among the various types of stars that populate the core of 
M15 we have identified: 



4.1 Blue straggler stars 

Blue straggler (BS) stars have bluer colours and are brighter 
than main sequence turnoff stars in globular clusters. The 
BS stars are frequently located in the densest regions 
of globular clusters where crowding makes difficult their 
identification. Previous work on M15 core stars has pro- 
duced BS candida tes selected on the basis o f CMDs in the 
the far-ultra v iolet iFerraro fc Parescd (1 1993h . I vs. U — V 
lYannv et all Jl994h or V vs. B - V llMoretti et all ()2008h 
using the data in IPiotto et all (|2002h ). Previously known 
blue stragglers are marked in figure fTTI with filled diamonds. 

In tabl e [3] we list the BS stars found in the p resent 
work and in I Yanny etal] i| 19941 ); iMoretti et alj ((2008) . Col- 
umn one gives the ID number in our catalogue, the following 
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Figure 12. Colour-magnitude diagram I vs. V — I for the M15 
core stars, /-band magnitudes are from FastCam (this work), V 
magnitudes are from HST (van der Marel et al. (2002)). Blue 
stragglers (BS), bright red giant branch (RGB), (red) horizontal 
branch (HB) and blue horizontal branch (BHB). Blue straggler 
stars given in previous references are marked with filled diamonds. 
Variable stars are taken from Dieball et al. (2007). 
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Figure 13. Location of blue straggler stars (BS) in the core of 
M15. The area showed is the field observed by FastCam. The co- 
ordinates (R.A., Dec.) for each star are measured in arcseconds 
with respect to star AC 211, marked with the + symbol in the fig- 
ure. The centre of the cluster given in van der Marel et al. (2002) 
is marked with a red cross. The SXPHE candidate is marked with 
a green circle. 



4.2 LMXBs 



three columns list the BS status according to each reference, 
and the last column gives some comments. The horizontal 
lines separate the previously known BS stars (upper part), 
new BS stars identified in this work (middle) and BS star 
candidates in previous works which are not confirmed as 
such here (bottom section). 

We find 10 BS candidate stars (location plotted in Fig- 
ure [13]), 7 were previously known and 3 are new candidates. 
One of the three new BS candidates is located at less than 
1 arcsec (r < 1") from the cluster centre. Three out of the 
seven previously known BS stars also have a distance from 
the cluster centre of less than 1.1 arcsec (r < l".l). 

We find in previous works 10 BS candidates which are 
not confirmed as such here. Only one of these 10 stars is 
listed in more than one reference (ID 145). For this object, 
we find a —1.41 mag difference betwe en the U magnitude 
by Yannv et al.l (| 19941 ) and the one by Ivan der Marel et al.l 
(2002). In the V-band the difference is —0.8 mag and com- 
paring the /-band in lYannv et all 1 19941 ) with ours the dif- 
ference is 0.41 mag. This ma y suggest variability, but the 
star is not listed as variable in lDieball et al] (|2007T ). 

The other not confirmed BS stars are listed as can- 
didates in only o ne previous work, the three found by 
lYannv et al] (|1994| ) are at less than 1 arcsec (r < 1") from 
the cluster centre where accurate photometry is more diffi- 
cult to perform due to crowding . The BS star candid ate ID 
369 is a variable star (ID V41 in lDieball et~all (|2007l )), clas- 
sified as CV, it is likely a SX Phoenicis (SXPHE) sta r This 
object is not listed in the catalogue of lYannv et al] (|1994| ) 
possibly because of its variability. 



Located in the inner region of M15 the V~ 15 star AC 211 
is optically one of the brightest known low mas X-ray binar 
(LMXB) systems. It was identified bv lAuriere et al.1 (|l98- 
as the possible optical count e rpart of the X-ray source 4U 
2127+119 and ICharles et all (|l98q) provided spectroscopic 
evidence for the classification as LMXB. The high optical to 
X-ray luminosity ratio suggests that a very luminous cen- 
tral X-ray source is hidden behin d the a ccretion disk. A de- 
tailed analysis bv lllovaiskv et al] (JT993) revealed an orbital 
period of 17.1 hr. This object is clearly identified in our Fast- 
Cam images (see the left panel of figure [T4")) as a source of 
/ = 16 (ID 102). We have searched for photometric variabil- 
ity grouping our series of data into 20 consecutive cubes (see 
the subsection on variability below). For each cube one final 
image was produced and photometry performed. The stan- 
dard deviation of the series of 20 photometric measurements 
is cr=0.05 mag, consistent with the statistical error. 

A seco nd luminous LMXB in th e core of M15 was dis- 
covered by I White fc Angelinil l|200ll) using the Chandra X- 
Ray Observatory. The Chandra observations resolved 4U 
2127+119 into two X-ray sources, the previously known 
AC 211 and a second source named as M15 X-2 (CXO 
J212958. 1+121002). This new X-ray source is 2.5 times 
brighter than AC 211 according to the Chandra counts 
rate and was a ssociated with a U = 18. 6 mag star in 
the data from [ Cuhathakurta et"afl (1 19961 ) (star 590 in 
Ide Marchi fc Parescd (Il994t) ) located at 3".3 from the M15 

report time-resolved FUV pho- 



Dieball et al 



centre. 

tometry for the optical counterpart and argue that the FUV 
emission is dominated by an irradiated accretion disk around 
the neutron star primary. These authors concluded that M15 
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Figure 14. FastCam /-band convolved images. Left panel: 
AC211 (M15 X-l), the yellow circle has a diameter of 0.5 arcsec 
and is centred in the AC211. Right panel: M15 X-2. The yellow 
circle is centred in the Chandra location of the X-ray detection 
and its diameter indicate the Chandra 0.5 arcsec error. The small 
white circle of diameter 0.14 arcsec marks the position of the FUV 
star associated to M15 X-2. 



Table 3. Blue straggler candidates. 



ID 


This work 


Yanny et al. 


Moretti et al. 


Comments 


263 


Yes 


Yes 


Yes 




276 


Yes 


Yes 


Yes 




285 


Yes 


Yes 


Yes 




320 


Yes 


Yes 


Yes 




239 


Yes 


Yes 


No 


r < 1" 


128 


Yes 


No 


Yes 


r < 1" 


369 


Yes 


No 


Yes 


r < I'M 


402 


Yes 


No 


No 


r < 1" 


156 


Yes 


No 


No 




337 


Yes 


No 


No 




145 


No 


Yes 


Yes 




103 


No 


Yes 


No 


r < 1" 


180 


No 


Yes 


No 


r < 1" 


220 


No 


Yes 


No 


r < 1" 


406 


No 


No 


Yes 




300 


No 


No 


Yes 




552 


No 


No 


Yes 




452 


No 


No 


Yes 




615 


No 


No 


Yes 




683 


No 


No 


Yes 





X-2 can be classified as an ultracompact X-ray binary. On 
the right panel of figure[l4]we show the location of the X-ray 
source in our final FastCam image and indicate the position 
of the FUV star. We set a 3cr upper limit to the I-band 
magnitude of this FUV source of I = 20.5. 



4.3 RR Lyrae, Cepheids and cataclysmic variable 
stars 

RR Lyrae stars (IDs 43, 48, 57, 68, and 75, in our catalogue) 
and Cepheid stars (I Ds 27, 44 and 46) previously reported by 
iDieball et all |2007l ') are marked with green circles in figure 
1111 We also detect the cataclysmic variable star (ID 369) 
which has colours suggestive of a SX Phoenicis (SXPHE) 
star a nd note that the othe r known cataclysmic variable star 
(V39, IDieball et all feOOTb ) with a magnitude of J ~ 21 is 
very marginally detected (2a level) in our image. This object 
seems to have a rather blue colour when compared with 



previous observations at shorter wavelengths (V = 18.65 
and U = 19.08). 



4.4 Variability: light curves of selected stars 

We have used our series of images to study the short 
timescale variability of stars of the previous types with / 
magnitude in the range 14.5-16 mag. For each set of 1000 
images (each obtained over a period of 30 s real time) we 
generated a lucky image. 200 such images were produced, 
each 10 consecutive images were then averaged in order to 
increase S/N and perform photometry on the stars of in- 
terest. The total time span for the resulting series of 20 
images is of 2 hours 43 min, each image corresponding to a 
time interval of 8.1 min. The series of photometric measure- 
ments obtained for a selection of the RR Lyrae, Cepheid, 
and LMXB stars in our image are shown in figure 1151 The 
majority of these stars show very little variability in the I- 
band. We note smooth trends in several of them (ID 43, 44, 
46) which are consistent with the known periods of variabil- 
ity. In some cases, we also detect changes among consecutive 
8.1 min images which are clearly beyond the measurement 
errors and may be indicative of intrinsic short time variabil- 
ity, for instance due to pulses. Remarkably, the LMXB star 
AC 211 did not show evidence for any significant variability 
on these time scales. 



5 CONCLUSIONS 

Observations of the M15 core with the 2.5-m Nordic Opti- 
cal Telescope using FastCam and the lucky imaging tech- 
nique provided /-band images with spatial resolution and 
sensitivity close to those obtained by HST in this region. 
We cross-match with the published HST/WFPC star cat- 
alogues and calibrate photometrically and astrometrically 
our image. The number of objects recovered in our final 
13" x 13" FastCam image is comparable to that reported 
for the HST images. A catalogue of 1181 stars is presented 
in this paper. Based on number counts the limiting mag- 
nitude of the catalogue is I « 19.5. Errors in magnitudes 
and positions are estimated from simulations. These simu- 
lations also indicate that crowding and spatial resolution, 
more than sensitivity, limit the completeness of the cata- 
logue to about one magnitude brighter. Using a CMD V 
vs. U — I we discuss the various stellar populations present 
in the M15 core. In particular, we identify a few new Blue 
Straggle star candidates and several new core stars which 
were not previously reported. We show that this imaging 
technique is particularly useful to investigate stellar popu- 
lations in the core of globular clusters where the presence 
of bright stars may limit the use of more conventional tech- 
niques. Lucky imaging observations of the core of Ml 5 and 
other globular clusters undertaken with a baseline of several 
years may provide proper motions of stars in the very inner 
region of the clusters (r < 1") with the precision required to 
constrain the properties of intermediate-mass black holes. 
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Figure 15. Mean-subtracted light curves (mag) for the variable sources previously identified as such in Dieball et al. (2007). There 
are 20 consecutive images which are represented in the x axis. Observations started at BJD-2454672.5215 and each image corresponds 
to a time interval of 8.1 min. The RR Lyrae stars arc IDs 48, 75, 68, 43 and 57 and they are identified as V9, V14, V23, V28 and V29 
respectively in Dieball et al. (2007). The Cepheid stars are IDs 27, 44 and 46 and they are identified as V10, V13 and V18 respectively 
in that reference. ID 102 is AC 211 and it is identified as V21 and ID 42 is identified as V24 in Dieball et al. (2007), where they found 
V24 was in the BHB zone but we find it is in the HB zone. Note that the scale of the upper curve is larger than the others. Errors are 
at 1-cr. 



ACKNOWLEDGMENTS 



We thank the Nordic Optical Telescope staff and the Instru- 
ment Maintenance team of the IAC for their support dur- 
ing the observations. Based on observations made with the 
Nordic Optical Telescope operated on the island of La Palma 
by the NOTSA in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astroffsica de Canarias. This 
research has been supported by by Project No. 15345/PI/10 
from the Fundacion Seneca and the Spanish Ministry of 
Economy and Competitiveness (MINECO) under the grant 
AYA2010-21308-C03-03. This research has made use of the 
VizieR catalogue access tool, CDS, Strasbourg, France. 
Based on observations made with the NASA/ESA Hub- 
ble Space Telescope, obtained from the data archive at the 
Space Telescope Science Institute. STScI is operated by the 
Association of Universities for Research in Astronomy, Inc. 
under NASA contract NAS 5-26555. 



REFERENCES 

Auriere, M., Le Fevre, O., Terzan, A, 1984, A& A, 138, 
415. 

Baldwin, J. E., Warner, P.J., Mackay, C. D, 2008, A& A, 
480, 589. 

Charles, P. A.; Jones, D. C; Naylor, T., 1986, Natur, 323, 
417. 

Davidge, T.J. & Courteau, S., 1999, AJ, 117, 1297. 
Debray, B., Llebaria, A., Dubout-Crillon, R., Petit, M., 

1994, A& A, 281, 613. 
de Marchi, C, Paresce, F., 1994, ApJ, 422, 597. 
Dieball, A., Knigge, C, Zurek, D. R., Shara, M. M., Long, 

K. S., Charles, P. A., Hannikainen, D, van Zyl, L., 2005, 

ApJ, 634, 105. 

Dieball, A., Knigge, C, Zurek, D. R., Shara, M. M., Long, 
K. S., Charles, P. A., Hannikainen, D., 2007, ApJ, 670, 
379. 

Femema, B., Rebolo, R., Perez-Prieto, J. A., Hildebrandt, 
S.R., Labadie, L., Peerez-Garrido, A., Beejar, V.J.S., 
Diaz-Sanchez, A., Villo, I., Oscoz, A., Lopez, R., 
Rodriguez, L.F., Piqueras, J., 2011, MNRAS, 413, 1524. 

Ferraro, F.R. k Paresce, F., 1993, AJ, 106, 154. 



© 2010 RAS, MNRAS 000, 000-000 



High-resolution optical imaging of the core of the globular cluster Ml 5 with FastCam 



Gebhardt, K., Pryor, O, O'Connell, R. D., Williams, T. 
B., Hesser, James E., 2000, AJ, 119, 1268. 

Guhathakurta, P., Yanny, B., Schneider, Donald P., Bah- 
call, J.N., 1996, AJ, 111, 267. 

Ilovaisky, S. A., Auriere, M., Koch-Miramond, L., Cheva- 
lier, O, Cordoni, J.-P., Crowe, R. A., 1993, A& A, 270, 
139. 

Johnson, H.L., 1966, ARA&A, 4, 193. 

Kulkarni, S. R., Goss, W. M., Wolszczan, A., Middleditch, 
J., 1990, ApJ, 363, L5. 

Labadie, L., Rebolo, R., Villo, I., Perez-Prieto, J. A., 
Perez-Garrido, A., Hildebrandt, S.R., Femema, B., Dfaz- 
Sanchez, A., Beejar-Sanchez, V., Oscoz, A., Lopez, R., 
Piqueras, J., Rodriguez, L.F., 2011, A& A, 526, 144. 

Law, N. M., Mackay, C. D., Baldwin, J. E., 2006, A& A, 
446, 739. 

Mason, B.D., Wycoff, G. L., Hartkopf, W. I., Douglass, G. 

G., Worley, C. E., 2001, AJ, 122, 3466. 
McNamara, B. J., Harrison, T. E., Anderson, J., 2003, ApJ, 

595, 187. 

Moretti, A., de Angeli, F., Piotto, G, 2008, A& A, 483, 183. 

Oscoz, A., Rebolo, R., Lopez, R., et al. 2008, in Presented 
at the Society of Photo-Optical Instrumentation Engi- 
neers (SPIE) Conference, Vol. 7014, p. 701447, Society 
of Photo-Optical Instrumentation Engineers (SPIE) Con- 
ference Series. 

Piotto, G.; King, I. R.; Djorgovski, S. G.; Sosin, C; Zoccali, 
M.; Saviane, I.; De Angeli, F.; Riello, M.; Recio-Blanco, 
A.; Rich, R. M.; Meylan, G.; Renzini, A., 2002, A& A, 
391, 945. 

van der Marel, R.P., Gerssen, J., Guhathakurta, P., Peter- 
son, R.C., Gebhardt, K., 2002, AJ, 124, 3255. 

White, N. E., Angelini, L., 2001, ApJ, 561, 101. 

Yanny, B., Guhathakurta, P., Bahcall, J.N., Schneider, 
D.P, 1994, AJ, 107, 1745. 



© 2010 RAS, MNRAS 000, 000-000 



